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Abstract 

We extend the Zee model, where tiny neutrino masses are generated at the one-loop level, to a 
supersymmetric model with R-parity conservation. It is found that the neutrino mass matrix can 
be consistent with the neutrino oscillation data thanks to the nonholomorphic Yukawa interaction 
generated via one-loop diagrams of sleptons. We find a parameter set of the model, where in 
addition to the neutrino oscillation data, experimental constraints from the lepton flavor violating 
decays of charged leptons and current LHC data are also satished. In the parameter set, an 
additional CP-even neutral Higgs boson other than the standard-model-like one, a CP-odd neutral 
Higgs boson, and two charged scalar bosons are light enough to be produced at the LHC and future 
lepton colliders. If the lightest charged scalar bosons are mainly composed of the SU(2)i-singlet 
scalar boson in the model, they would decay into en and go with 50 % of a branching ratio for each. 
In such a case, the relation among the masses of the charged scalar bosons and the CP-odd Higgs 
in the minimal supersymmetric standard model approximately holds with a radiative correction. 
Our model can be tested by measuring the specific decay patterns of charged scalar bosons and 
the discriminative mass spectrum of additional scalar bosons. 
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I. INTRODUCTION 


Neutrino oscillation data [l-9| have indicated the existence of tiny masses of neutrinos, 
which are absent in the standard model (SM) of particle physics. If the tiny neutrino masses 
are generated by a new physics at a very high energy scale (e.g., the seesaw mechanism (l0|), 
such a new physics (heavy new particles) is not directly accessible by experiments. In 


contrast, in scenarios based on radiative generation of tiny neutrino masses 
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16|, the 


smallness of neutrino masses is deduced by the quantum effect without introducing very 
heavy new particles. Therefore, such scenarios can consequently be tested by current and 
future collider experiments. 

The model proposed by A. Zee [n| is the first model of radiative generation of neutrino 
masses, where an SU(2)L-singlet scalar held with a hypercharge Y = 1 and the second 
SU(2)i-doublet Higgs held with H = 1/2 are introduced to construct the one-loop diagram 
for the neutrino mass. Studies on the phenomenology in the Zee model can be found in 
In the minimal version of the model (the so-called Zee-Wolfenstein model ji^h 


Refs. 


17 


lepton havor violating (LFV) Yukawa couplings with the second Higgs doublet are forbidden 
at the tree level. However, such a model has already been excluded by current neutrino 
oscillation data (See e.g.. Ref. [2^). In order to reproduce the neutrino data, lepton havor 


violating interactions are necessary 


22 


24 


m. 


In the Zee model, there are several problems to be solved. Although the Zee model 
with LFV couplings is phenomenologically acceptable, they should be well controlled by 
some principle in order to suppress the dangerous havor changing neutral current (FCNC) 
processes^ The model is also confronted by the quadratic divergence problem like the 
SM. In addition, there is no dark matter (DM) candidate in the Zee model. If we consider a 
supersymmetric (SUSY) extension of the Zee model, we may be able to solve these problems 
simultaneously. The quadratic divergence is automatically cancelled by the loop contribution 
of SUSY partner particles. If the R-parity is conserved, the lightest SUSY particle becomes 
stable and it can be DM. Moreover, the LFV Yukawa couplings can naturally be induced. 


^ It would be also required to introduce the third SU(2)i-doublet scalar field which is only for the quark 
Yukawa interactions without FCNC in the quark sector. 

^ An idea to control such FCNC well so that stringent constraints on ^ —>■ eee and ^ ^ are automatically 
satisfied is introducing the A 4 symmetry to the Zee model 0. 
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The previous study of the SUSY extension of the Zee model is found in Ref. |20|, where 
the R-parity violation is introduced to the Minimal SUSY SM (MSSM). In this model, 
right-handed sleptons play the role of the charged singlet scalar in the Zee model. Since 
the sleptons carry lepton flavors in contrast with the singlet scalar in the Zee model, the 
flavor structure of the generated neutrino mass matrix becomes different from the one in the 


Zee-Wolfenstein mode 
can be found in Refs. 


. SUSY models for the other scenarios of radiative neutrino masses 


27 




In this paper, we propose a SUSY extension of the Zee model with the R-parity conser¬ 
vation (the SUSY Zee model). The stability of the dark matter candidate is guaranteed. 
This is the simplest alternative example of the MSSM with the right-handed neutrino su¬ 
perfields (the SUSY seesaw model). In the SUSY Zee model, the MSSM is extended by 
introducing a pair of SU(2)L-singlet superfields with hypercharge Y = 1 and —1, which 
also carry lepton numbers. The Higgs sector of the MSSM is the type II two Higgs doublet 
model [Mi at the tree level. The extra Higgs doublet can play the role of the second Higgs 
doublet of the Zee model. In SUSY models, nonholomor phi c Yukawa interactions are gen¬ 
erally induced by the one-loop effect of SUSY particles 32H37i . This mechanism may be 
utilized for generating the flavor violating interaction, which is required for the Zee model 
to satisfy the neutrino data. The structure of the LFV Yukawa matrix is determined by the 
flavor structure of the soft SUSY breaking slepton mass matrices. Such radiatively induced 
coupling constants are expected to be much smaller than the other Yukawa coupling con¬ 
stants. In this model, we study the neutrino mass matrix, and we find a benchmark point 
for model parameters, which satisfies the required structure of the neutrino mass matrix and 
the constraints from experimental searches for the LFV decays of charged leptons. On the 
benchmark point, phenomenological consequences of our model are discussed, and testabil¬ 
ity of our model at the LHC and future lepton colliders such as the International Linear 
Collider (ILC) is mentioned. 

This paper is organized as follows. In Sec. [TTl our model is defined. Section HTTI is devoted 
to showing how the neutrino masses are generated at the one-loop level. In the section, 
a benchmark set of model parameters which satisfy neutrino oscillation data are given. 
Phenomenology in the SUSY Zee model with the benchmark set is discussed in Sec. IIVI 
Conclusions are given in Sec. |Vl 
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TABLE I: Superfields of the SUSY Zee model. The baryon number is zero for all particles in this 
table. 

II. THE MODEL 

Superfields of the SUSY Zee model are partially listed in Table HI The transformation 
property under the R-parity is given by (—where B (L) is the baryon (lepton) 
number and s denotes the spin. Component fields with ’’tilde” are odd under the R-parity. 
The relevant part of the superpotential is constructed as 

W =y^tLj{-ia2)^d+ ^2 > ( 1 ) 

where \ (* = 1"3) are the Pauli matrices. The relevant part of the 

soft-SUSY breaking terms is given by 

-Fh.c.} 

- 1 ^ 12^2 ^2 - {Be, + h.c.} 

- {m~)u'L\Le - 

~ ^y£{.^E)ee> Lj,{ia 2 )^d + ^-c.^ — L£{ia 2 )L^i uf -|-h.c.| 

-|Ci$I$da;+ + C2$i<f)„a;2- + h.c.} , (2) 
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where = —A^^. Notice that the term of Ci gives the important interaction reqnired in 
the non-SUSY Zee model as the source of the lepton number violation by two units.^ 

In order to generate the neutrino mass matrix at the one-loop level, the following Yukawa 
interactions are used: 


Yukawa — He Lji{—icr2)^d + (h2)a' 

+ Aa^)££' LJ {ia2)L£i uf + {Y2a)u' LJ {ia2)Lii , 


(3) 


where the first and the third terms are obtained from the superpotential W in eq. ([T]) . The 
second and the fourth terms are generated at the one-loop level |32l-l37l|. In eq. (jS]), Ya 2 can 
be ignored, when is very heavy. We here consider such a case. The Yukawa matrix Y 2 is 
generated through the slepton mixing as follows 32l-l37l|: 


{Y2)te.’ — yA [s-AAii' + \i 




0!2 


2 / 3 (Mi , , m~ ) + / 3 (Mi , m~ ) - 2 l 3 (M ^, m~ ) 




^-y,M 2 

|■/i^M2(Am|) 




(4) 


(5) 


02 


££' 


££' 


'■L '-R '-L 

,2 


2 / 4 (Ml, rrij , m~ , ) -h / 4 (M^, ) 

2 h{M^,mj ,m| ,m| ) + 2 h{Ml, ,m| ) 


h{M ^, /i^,, mj ,mj,) + 2 h{M ^, /i^, m~ , m~ ) 




( 6 ) 


where Mi (M 2 ) is the soft SUSY breaking mass of U(l)y (SU(2)i) gauginos. The matrices 
Am- and Am- denote off-diagonal parts of {mY)££' and respectively. Thus, ( 62 )m = 

0. Loop functions l 3 {x,y,z) and l 4 {x,y, z,w) are defined as 


^ _ xy\n{x/y) + yz\ia.{y/z) + zx\n{z/x) 


h{x,y,z,w) = 


{x - y){y - z){z - x) 

X In X 1 / In 1 / 


(7) 


{y - x){z - x){w - x) {x - y){z - y){w - y) 


^ Such three-scalar terms with both fields and conjugate fields are the so-called C-terms [3^, which are 
usually ignored because most SUSY breaking scenarios do not generate the C-terms. However, the C- 
terms can be generated in some SUSY breaking scenarios, e.g., in an intersecting D-brane model with the 
flux compactification [39|. Some detailed discussion about the SUSY breaking effects for the radiative 
neutrino mass can be found, e.g., in Ref. 4Cll |. 
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zlnz wlnw 

(x — z){y — z){w — z) (x — w){y — w){z — w) 

Although all terms in eq. (E]) (in eq. ([5])) are proportional to /i^, the first and the third 
terms in eq. ([H]) (the first term in eq. dS])) do not contain Higgsinos in the loop. Therefore, 
sizable (and ei also) can be generated by these terms if we take much larger /i<j, than the 
other mass scales (e.g., Mi, rn~). Then, ei and 62 are almost independent of the value of 
M 2 . Yukawa interactions in eq. ([3]) can be rewritten as 

Yukawa = -- LJ H - - Xi£i + 0^a)££' LJ {ia 2 )Lii ujf , (9) 

V V 

where + nj = (246 GeV)^, = \/2(0°), and Vd = \/2(0j)). The matrix Ya is an 

arbitrary antisymmetric matrix. Two Higgs doublet fields $o and are defined as 



/ Cfj -sp 

\sp Cp 


(-icT2)<hl 


( 10 ) 


where sp = sin/3 and cp = cos/3 for tan/3 = Vu/vd- Interactions in eq. (jO]) are expressed in 
terms of mass eigenstates of charged leptons^, which are given by diagonalizing the mass 
matrix cpy^ + spY 2 . By ke^ng ei and 62 up to linear terms (keeping ei tan/3 for all order), 
the matrix X is given by 321 

1 + tan^ /3 


^£e' — ~ tan /3 5££i + 


(1 + tan/3(ei)£)2 


{^l)£Si£> + {e2)l£i 


( 11 ) 


The off-diagonal elements of X provide FCNC, and they are important to obtain the appro¬ 
priate structure of the neutrino mass matrix. 

Since Nambu-Goldstone modes are contained in mass eigenstates of the three charged 
bosons are given by linear combinations of uf, and a charged component (pQ of $o- The 
matrix of squared masses of charged bosons is given in a basis of (0 q , 0 ;/*“, cu^) by 



-h m\ 

Civ 

v/2 

C2V \ 

^/2 

M^, = 

Civ 

h2 

{Mjj+)22 

Bu]y^ 


C 2 V 

\ %/2 


(TfH+)33 y 

(M^+)22 = - 

-m^ tan^ dw cos(2/3) -|- 


{MhPss = m/p tan^ cos(2/3) + m/, 2 + ;i/, 


( 12 ) 

( 13 ) 

( 14 ) 


^ Charged leptons in eq. (|3]) are not mass eigenstates although we used the same notation i. 


6 













u 


1 / 


Cl 

M, 


\ 'TO 


Ya 




£" 


R 


^MfX 


^I'L 


FIG. 1: The one-loop diagram for light Majorana neutrino masses in the SUSY Zee model. 


where tua is the mass of the CP-odd Higgs boson A^. The matrix Mjj+ is diagonalized as 
= [/jif+diag(m^+,with a nnitary matrix Uh+- We here assnme that 
(M ^+)33 is much larger than the other elements for simplicity, so that mixing effects via C 2 
and can be ignored. The mass eigenvalues of charged scalar bosons are given by 



The mixing matrix is given by 

(. 

Uh+ = 


COS0_| 

sin 6^4 
0 


- sin 9a 0 
cos 9 a 0 
0 1 


2„,2 


sin^ 29+ = 


2Cfv 




+ m\- m^4 




rrij. — m. 

Hr, 




(18) 


(19) 


For 9a — 0, the charged scalar boson is the singlet-like one ~ o;]^) while is 
almost the same as the charged Higgs boson of the MSSM. 


III. NEUTRINO MASS AND BENCHMARK SCENARIO 

The neutrino mass matrix {m^)ui in the flavor basis, (l/ 2 )(m^)£^/(i 2 £j^)'=z/^,j^ -|- h.c., is 
generated by the diagram in Fig. [T]and its transpose diagram. Keeping the leading terms 
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of the charged lepton masses, the neutrino mass matrix is expressed as 


- 

f^loop 

V 


YaM}X + (YaM}X)'‘ 


( 20 ) 


where Mi = diag(me, m^, 171 ^-) and Cioop is a flavor-independent factor calculated by the loop 
integration. The explicit form of Cioop is given by 

,2 

( 21 ) 


m. 


^ sin 20 + 

^loop —-^ 


327r2 


m 




The case with X = 1 corresponds to the Zee-Wolfenstein model. Since we assume rather 
heavy SUSY particles in order to satisfy constraints from LFV decays of charged leptons, 
we can ignore a contribution from the loop.^ 

A mass matrix for Majorana neutrinos in the flavor basis can be diagonalized by 

the Maki-Nakagawa-Sakata (MNS) matrix ^ Umns as 


= ClMNsd‘‘is(™ie”“. ™2, m3e'“”)flit,Ns. 


( 22 ) 


where a ^2 and a ^2 are Majorana phases The MNS matrix can be parametrized as 


Umns — 


10 0 


Ci3 0 Si 3 e 


C 12 S 12 0 

0 C 23 S 23 


0 1 0 


-S 12 C 12 0 

1 0 —S 23 C 23 , 


1 ~Si 3 e*^ 0 Ci 3 j 


0 

0 


(23) 


where Cj,- = cos 0^2 and Sj,- = sin0j,'. The measurement of the u disappearance at the T2K 


experiment jo] shows 


sin^ 023 = 0.514l°;°“, Aml 2 = (2.51 ± 0.10) x 10"^ eV^ (24) 

for the normal mass ordering (mi < m 3 ) and 

sin^ 023 = 0.51ll|(;J(|, Am^g = (2.48 ± 0.10) x 10"^ eV^ (25) 

for the inverted mass ordering (m 3 < mi). A combined analysis of the solar neutrino 
measurements and the KamLAND data results in 


tan^ 012 = 0.427lK;[(^I, Am^i = 7.46;°;^° x 10 "^ eVl (26) 


® The interaction wj*" in Fig.[T]is replaced with a nonholomorphic Yukawa interaction uj^, which 
is generated via a one-loop diagram with $ 0 , ^nd the bino. 
















The reactor measurement at the DayaBay experiment 8 | gives 

sin^ 2013 = 0.084 ±0.005. (27) 

Let us define (YA)efj, = and keep only terms of in eq. (I20|) . We then 

obtain = 0 , which can be satisfied for the inverted mass ordering with the following 

values: 

sin^ 023 = 0.511, sin^ 20i3 = 0.09, tan^ 0i2 = 0.427, (28) 

ai 2 = 032 = tt, 6 = 0, (29) 

= 2.48 X 10"^ eV^ Amj^ = 7.46 x lO"'^ eV^ m 3 = 5.05 x lO^^ eV. (30) 

For these values, the structure of the neutrino mass matrix is given by 


^-2.86 4.23 -4.81^ 


m„ = 


X 10 ^eV. 


(31) 




4.23 -2.14 -3.91 
y-4.81 -3.91 0 

When we search for a set of model parameters which gives the structure of m^ in eq. (I3TD . 
constraints from LFV decays of charged leptons have to be taken into account. In Table UTl 
we summarize current data from various LFV experiments. 

By assuming = 0, = Xt-t, and X^r = the neutrino mass matrix in eq. (l2n|) 

is determined by five combinations of model parameters: Xr^jX^-r., Xt^/X^- r, (XA)eT/(XA)fiT, 
(^A)eM/(^A)/.r, and Cioopm‘^{YA)f,rXrr‘ They are constrained as 

(Td)er ^re 






(n)eM X,, 


(Y. 


(r;) 


A) Xtt 

X. 


± 




(lA)Mr {YA)>.r X. 


V 2 , 


iYA)er 
iYA),r 

Tfi ^re 

Xr, 



(mu] 

ler 

{mA 

'flT 

{rnAj 


{rriAj 

'flT 

{rriu 





2 {my) 


= 0.366, 
= 1.23, 

= -1.08, 
= 0.273, 


flT 


—Cloov^AYA)^lTXrr = = -3.91 X 10 ^eV. 


(32) 

(33) 

(34) 

(35) 

(36) 


.6. 


These constraints are satisfied with the following benchmark set for model parameters” 

Mq = m~ = m~ = Ml = 10 TeV, M 2 = 3TeV, (37) 


® If the contribution from {YA)e^i to rrii, is naively ignored, larger off-diagonal elements of X are necessary. 
We do not take this option in order to suppress LFV decays of charged leptons as much as possible. 


9 






















Current bound 

Euture sensitivity 

BR(/r —>• ey) 

< 5.7 X 10-13 (MEG) [4^ 

6 X 10-^^ (MEG upgrade) [44] 

BR(// eee) 

< 1.0 X 10-12 (SINDRUM) [4^ 

~ 10-3® (Mu3e) [4^ 

BR(r ey) 

< 1.2 X 10-7 (Belle) [47] 

< 3.3 X 10-s (Babar) [4^ 


BR(r —>■ /ry) 

< 4.5 X 10-8 (Belle) [47] 

< 4.4 X 10-8 (Babar) [48] 

5 X 10-3* (Belle II) [4^ 

BR(r —)• eee) 

< 2.7 X 10-8 (Belle) [^ 

< 2.9 X 10-8 (Babar) [^ 


BR(r — )■ ee^) 

< 1.8 X 10-8 (Belle) [^ 

< 2.2 X 10-8 (Babar) [^ 


BR(r —)• e/i/i) 

< 1.7 X 10-8 (Belle) [^ 

< 2.6 X 10-8 (Babar) [^ 


BR(r ^ /lee) 

< 1.5 X 10-8 (Belle) [^ 

< 1.8 X 10-8 (Babar) [^ 


BR(r —)• /le/x) 

< 2.7 X 10-8 (Belle) [^ 

< 3.2 X 10-8 (Babar) [^ 


BR(r —)• /I/x/i) 

< 2.1 X 10-8 (Belle) [^ 

< 3.3 X 10-8 (Babar) [^ 

< 4.6 X 10-8 (LHCb) [^ 

1 X 10-3* (Belle II) [4^ 

BR(r ^ fj-T]) 

< 6.5 X 10-8 (Belle) [^ 



TABLE 11: Constraints and future sensitivities for LEV decays of charged leptons at 90 % confi¬ 


dence level (CL). 

= -(0.8 Mo)^ (38) 

(Am|),e = (Amf),e = -(0.708 M^)\ (39) 

(Am|)^e = (Am|)^e = 0, (40) 

/i$ = 762 Mq, tan/3 = 2, (41) 

mA = 380 GeV, rn„+ = 350 GeV, sin^6*+ = 10“^, (42) 

-^1 
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iYA)^r = -1.31 X 10-^ {YA)er = -2.24 X 10 


-4 


m 


(lA)eM = = 6.50 X 10 


-1 


m 


( 43 ) 

( 44 ) 


where ~ m\ + = (388 GeV)^ at the tree level. For the mixing angle a of CP- 

even neutral Higgs bosons, we obtain cos{(3 — a) = —0.027 from tan(2a)/tan(2/3) = -|- 

771 “^^)/{m\ — m\). Notice that a low tan/? value and large scales of and soft breaking 
parameters in the benchmark set are chosen in order to satisfy constraints from neutrino 
oscillation data and searches for the LFV decays of charged leptons. Nevertheless, even 
taking such a low tan/? value, it can also be compatible with mh = 125 GeV if we take 
M| = > (0(10) TeV)^ [5^, where m~ and m~ are masses of two stops, and such a 

large value of Ms is consistent with large breaking scales in the benchmark set. 

The value = 380 GeV satishes the constraint rriA ^ 350 GeV for tan/3 = 2 |55| 
which comes from the —)■ Zh^ search at the CMS experiment with the 19.7fb~^ data 
at y/s = 8 TeV 5g . The values at the benchmark point for rriA and tan /3 also satisfy the 
constraint at the ATLAS experiment with the 20.3 fb~^ data at y/s = 8 TeV 57|. The value 
m + = 350 GeV for ~ uf is consistent with a constraint on m~for a massless neutralino, 

Ul i 

which is 250 GeV obtained at the ATLAS experiment for 20.3 fb~^ of the integrated 

luminosity with y/s = 8 TeV j^. The CMS experiment gives > 200 GeV for 19.5 fb“^ 
of the integrated luminosity with y/s = 8 TeV 591 ]. 

The Yukawa matrix for <Fq in eq. ([9]) with the benchmark set is calculated as 


( 


X' = —M,X = 


-1.63 X 10-5 
0 

-1.22 X 10-2 




0 -3.49 X 10 

-3.39 X 10-3 9.25 x 10-^ 
-1.55 X 10-2 -5.69 X 10-2 


(45) 


/ 


This matrix is the source of FCNC. 


IV. PHENOMENOLOGY 

A. Lepton Flavor Violating Decays of Charged Leptons 

For ^ £'7, contributions from one-loop diagrams with charged scalars Hf are negligible 

because (V 4 )«//(V 4 )^//£/ are small enough. Since slepton masses are 0(10) TeV, one-loop 
diagrams involving sleptons have only negligible contributions to £ —)■ £'7 although 
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and {Am‘~)rii have a similar size to m| (also for m~). The condition in eq. fHOj) forbids not 
only the one-loop contribution of sleptons to ey but also the Barr-Zee type two-loop 


contributions 


60 


6l| with X'^^. The dominant contribution to /r —?■ ey comes from a one- 
loop diagram involving r with which results in BR(/i ey) = 1.5 x 10“^^ with 

niH ~ niA- This value satishes the current constraint BR(/i —)> ey) < 5.7 x 10“^^ (90% CL) 
at the MEG experiment and it could be observed at a planned MEG experiment 


upgrade 


44| where a sensitivity for BR(/i —)■ ey) ~ 6 x 10 is expected. If we take 


a larger value of tan/3 (e.g., tan/3 = 3), X'^^ is enhanced, and then it is required that 
(and also X/.g in fact) becomes larger due to the condition in eq. fl5^ : thus, a low tan/3 
value is required to satisfy the constraint on BR(/i ey) with = 380 GeV which is 
experimentally accessible. On the other hand, Barr-Zee diagrams involving the top quark 
in a loop give dominant contributions to r —)■ £y because {Am‘~)re, {Am~)re, {Am‘~)r^, 
and {Am~)r^ are not zero at our benchmark set. By using formulae in e.g. Ref. 6l| with 
BR(r — eUrVe) = 0.17, we have BR(r ey) = 2.7 x 10“® and BR(r —?■ /ly) = 4.4 x 10“® 
which satisfy BR(r ey) < 3.3 x 10“® (90 % CL) and BR(r —)■ /ly) < 4.4 x 10“® (90 % CL) 


obtained at the Babar experiment 48|]. These values are comparable to expected sensitivity 


BR(r i'y) ~ 10 ® at the Belle II experiment ^|. 


The Yukawa matrix X' can cause /i eee and r at the tree level. Although 

there is the stringent experimental constraint BR(/i eee) < 1.0 x 10“^^ (90% CL) at 
the SINDRUM experiment 0 . our benchmark set trivially satishes it because X'^^ = 0 in 
eq. (H5il gives BR(/i —)■ eee) = 0 at the tree level. For r the tree level contributions 

at the benchmark set result in BR(r — )■ e/x/i) = BR(r Jlee) = 0, BR(r eee) ~ 
BR(r —ee/i) ~ 10“^®, and BR(r — )■ fleii) ~ BR(r —)■ /Z/r/i) ~ 10“^^ where experimental 


bounds are BR(r —)■ ££'£") < 10 ® (90% CL) 50l-l52|. The constraint BR(r ^ ^rj) < 


6.5 X 10“® (90% Cld [53| is also satished even if we use the relation BR(r — ?■ /ly) ~ 
8.4 X BR(r —/J/U/i) 62|. 

The LEV coupling X'^^ = —1.55 x 10“^ in Eq. fjTHll is comparable to \/2mr/v = 10 
However, the branching ratio BR(/i —)■ /rr) is suppressed to about lO”"^ by cos^(a—/3) ~ 10 


-2 

-3 


Therefore, the benchmark point satishes the current upper limit on BR(h fir) 63|, l64j . 


although the 2.4 a excess which is currently reported by the CMS 63(] is not explained by 
our benchmark point. 
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FIG. 2: Decay branching ratios of (left) and H 2 (right) with respect to sin^0_|_. The red 
thick and the blue thin solid lines show BR(fl^ ev) (= —)• i^v)) and BK{H~ 

respectively. The magenta dashed line is for BK{H^ —)• tb). 


TV) 


B. Dark Matter 


Large values of and Mi are preferred in order to obtain sizable off-diagonal elements of 
X (namely, 62 ) which are required for the appropriate structure of the neutrino mass matrix. 
On the other hand, the value of M 2 is not required to be very large. Therefore, there is a 
possibility of the wino dark matter in the SUSY Zee model. In our benchmark set, we take 


M 2 = 3 TeV for which the relic abundance of dark matter can be explained 6 ^ . The spin- 
independent cross section of the wino scattering on a proton is evaluated as 10 ^’^cm^ (See 
e.g.. Ref. 6 ^), which is greater than the neutrino background 671. 


C. Phenomenology of Charged Scalar Bosons 

Decay branching ratios of Hf and H 2 are shown in Fig. [2] as a function of sin^ 6 *_|_ where 
Hf = uf if sin^ 6 ^+ = 0. The red thick and the blue thin solid lines show BR(if^ 
ev) (= BR(iL^ fiv)) and BR(iL^ — )■ tv), respectively. The magenta dashed line is for 
BR(Tr^ tb). Except for sin^ 0+ and Y^, parameters are set to the benchmark point. 
Values of elements of Ya depend on sin^ through a condition eq. fl36l) . For sin^ 9^ > 10“^, 
FTf is the doublet-like Higgs boson for which the decay —)■ tb is the dominant channel. 
In this case, since the dominant decay channel of H 2 is also tb (for any sin^ 9^), existence of 
an SU(2)-singlet component is hidden. Such doublet-like charged Higgs bosons (the same as 


13 

















the one in the Type-II two Higgs doublet model) with m ± = 350 GeV and m ± ~ 388 GeV 

^1 ^2 


for tan P = 2 may be observed at the LHG with y/s = 14 TeV and 300 fb~^ of the inte gra ted 
luminosity via the production process gb —)■ tH^ followed by the decay —)■ tb 681, iGOj. 

On the other hand, dominantly decays into leptons via Ya for sin^ 9^ < 10“^. The 
hierarchical structure of Ya in eqs. fjTHl) and fj441) gives a characteristic prediction 


BR{H^ ev) : BR{H^ fxu) : BR{H^ ^ rz/) ~ 1 ; 1 : 0, 


± 


(46) 


where neutrinos in the hnal states are summed because experiments are not sensitive to 
their flavors. If the coincidence BR{H^ eu) = BR{H^ —)■ /iz/) is observed, it would be 
regarded as a nonaccidental one but a natural consequence of the SU(2)/^-singlet charged 
scalar with a (YA)e/i-dominated Yukawa matrix. Nonobservation of the signal of tu 

would also suggest that Hf does not come from an SU(2)i-doublet Higgs which has a 
vacuum expectation value. For the singlet-like charged scalar, a region of < 430 GeV 


can be probed at the LHG with yfs = 14 TeV and 100 fb ^ of the integrated luminosity [70|. 


At the ILG with y/s = 1 TeV, the cross section is about 10 fb [2l| which would be enough to 

observe Hf. The singlet-like Hf would be distinguished from right-handed sleptons if the 

lightest neutralino is sufficiently heavy such that it does not mimic a neutrino. 

The behavior of m? + at the tree level with respect to sin^6'+ is shown in Fig. |3l In 
-^2 

the region of small sin^6'+, the relation rrp + = m\ + rnpy approximately holds with an 

nn 

appropriate radiative correction 7]|. On the other hand, the relation between and the 
mass of the charged Higgs boson in the MSSM at the one-loop level is expressed as 


m^+ =m\ + m%r + BAA{m‘^A) “ ^H+H-{mA + ^w) + 


(47) 


where formulae of self-energies B-AAig^), and Bww{.g^) can be found in Ref. [72|. 


This relation at the one-loop level also approximately holds for m? . in the SUSY Zee model 


for a small value of sin0_|_, where is almost the same as the charged Higgs boson i 


m 


the MSSM. The benchmark set gives ~ 0.1, where the is defined as = 

H. ^ ^ 






±). If the mass relation is experimentally confirmed in addition to eq. fl46ll 
^2 

which is a consequence of rriu in eq. (l20l) . the SUSY Zee model can be highly supported. 

The detection prospect of with tua = 380 GeV seems marginal for -A- tt at the LHG 
with yfs = 14 TeV 55, 68, Tsl but sufficient for A^ -A- Zh 55|. Discovery of the can be 


a but sufficient for A^ ^ Zh 

. The cross section for a process e+e“ -A- Z* ^ H^A^ -A ttbb 
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FIG. 3: The mass of at the tree level with respect to sin^ 0+ 
= 350 GeV. 


for 


380 GeV and 


for the at the ILC with 
to detect the signal. 


= 1 TeV is greater than 0.1 fb 



which would be sufficient 


V. CONCLUSION 

We have extended the Zee model to a SUSY model with the conserved R-parity. The 
MSSM has been extended by introducing SU(2)i-singlet complex superfields Qf and In 
order to generate three nonzero neutrino masses, the extension gives the simplest alternative 
to the SUSY seesaw model where three right-handed neutrino superhelds are introduced. 
Tiny neutrino masses are obtained at the one-loop level. We have shown that the mass 
matrix can be consistent with the current neutrino oscillation data thanks to the nonholo- 
morphic Yukawa interaction, which is dominantly generated by one-loop diagrams involving 
sleptons and the bino. We have obtained a benchmark point which satishes not only neutrino 
oscillation data but also constraints from LFV decays of charged leptons and the current 
LHC results. The parameter set is also consistent with mu = 125 GeV. The dark matter is 
stabilized thanks to the R-parity conservation in the SUSY Zee model, and we have found 
that the wino can be dark matter. 

While SUSY particles are rather heavy at the benchmark point, additional scalar 
bosons (iL°, Hf, and Hf) are light enough to be discovered at the LHC and future 
lepton colliders. If the lightest charged scalar bosons Hf are almost singlet-like due to a 
small mixing, their decay branching ratios have been predicted as BR(iLj’" —)■ ef) : BR(iLj’" —)■ 
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fiu) : —)■ Tu) = 1:1:0. On the other hand, the heavier charged scalar bosons Hf 

decay into tb. For such a small mixing case, a relation between m ± and in the MSSM 
remains because m ± is almost the same as the charged Higgs boson mass in the MSSM. 
Therefore, our model can be tested by measuring the specihc decay patterns of and 
and the discriminative mass spectrum of additional scalar bosons. 
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